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ABSTRACT 

Recent observations of high ionization rates of molecular hydrogen in diffuse inter- 
stellar clouds point to a distinct low-energy cosmic-ray component. Supposing that this 
component is made of nuclei, two models for the origin of such particles are explored and 
low-energy cosmic-ray spectra are calculated which, added to the standard cosmic ray 
spectra, produce the observed ionization rates. The clearest evidence of the presence of 
such low-energy nuclei between a few MeV per nucleon and several hundred MeV per 
nucleon in the interstellar medium would be a detection of nuclear 7-ray line emission 
in the range ~ 0.1 - 10 MeV, which is strongly produced in their collisions with the 
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interstellar gas and dust. Using a recent 7-ray cross section compilation for nuclear col- 
lisions, 7-ray line emission spectra are calculated alongside with the high-energy 7-ray 
emission due to ir^ decay, the latter providing normalization of the absolute fluxes by 
comparison with Fermi-LAT observations of the diffuse emission above = 0.1 GeV. 
Our predicted fluxes of strong nuclear 7-ray lines from the inner Galaxy are well be- 
low the detection sensitivies of INTEGRAL, but a detection, especially of the 4.4-MeV 
line, seems possible with new-generation 7-ray telescopes based on available technol- 
ogy. We predict also strong 7-ray continuum emission in the 1-8 MeV range, which in a 
large part of our model space for low-energy cosmic rays exceeds considerably estimated 
instrument sensitivities of future telescopes. 

Subject headings: gamma raysdSM - cosmic rays 



Introduction 



The component below about one GeV per nucleon of the interstellar cosmic-ray spectrum 
is probably the least well known. Direct observations are not conclusive because cosmic rays 
with these energies are effectively deflected by the action of the outstreaming solar wind such 
that the energy spectra measured inside the heliosphere are strongly altered with respect to the 
interstellar one. It results in an uncertainty of the particle flux which is increasing towards lower 
energies and below a few hundred MeV per nucleon the cosmic-ray composition and spectrum 
are essentially unknown. Indirect evidence for an important low-energy component comes from 
observations of the interstellar ionization rate. In particular, observations of the molecule in 
diffuse interstellar clouds point to a cosmic-ray induced ionization rate which exceeds by more than 
an order of magnitude the one which is calculated for the stan dard cosmic-ray flux thought to b e 
produced by diffuse shock acceleration in supernova remnants (jlndriolo. Fields. &: McCall I (|2009l ): 
Indriolo &: McCall I (|2012l )). However, neither the nature - electrons or nuclei - nor the spectrum of 
the low-energy component can be deduced from these observations. An independent argument for 
the existence of a significant component of low-energy cosmic-ray nucle i comes from the observed 
quasi-linear increase of Be with metallicity ([Tatischeff &: Kiener I (|201ll )). 



The most compelling indirect observation of this component would be a detection of nuclear 
7-ray lines from the galactic disk. In fact, the strongest lines from the most abundant nuclei are pref- 
erentially produced in the cosmic-ray energy range below a few hundred MeV in collisions of protons 
and a-particles with interstellar matter. They are expected to be the same as frequently observed 
from strong so lar flares, i.e. lines from the deexcitat i on of th e first few levels in ^^C, ^^O , ^^Ne, ^^Mg, 



28Si and ^^Fe (jRamatv. Kozlovskv. Lingenfelter I (ll979|)). Ilndriolo. Fields, fc McCaU I (|2009|) con- 
structed cosmic-ray spectra which produce the observed ionization rate in diffuse clouds by adding 
a low-energy component to the standard cosmic-ray spectrum. With that they estimated fluxes 
of the 4.4-MeV line of ^^C and the 6.1-MeV line of ^^O from the central radian of the Galaxy 
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close to the sensitivity limits of INTEGRAL. Besides strong narrow lines the nuclear line emission 
is composed of broad lines produced by the cosmic-ray heavy-ion component and of thousands of 
weaker lines which together form a quasi-continuum in the ~ 0.1 - 10 MeV range. 

In this paper, we calc ulate this total nuclear 7-ra,y line emission from the inner Galaxy by 



making use of the wo rk of llndriolo. Fields 



McGah 



(200' 



n ew de velopments in nuclear re- 



action cross sections (jMurphv et al. iBenhabiles et al. I (120111)') and rece r it ob se rvations 



of cosmic-ray induced high-energ y 7-ray emission with Fermi-LAT (jAbdo et al. I (j2009l ): IStrong 



(j201ll ): lAckermann et al. I (j2012l )). Finally, the possibilities for a detection of this emission with 



existing and eventual future space-borne 7-ray instruments are discussed. 



2. Galactic cosmic-ray spectra and composition 



Following the method of llndriolo. Fields. &: McCall I (|2009l ). the interstellar cosmic-ray (CR) 
spectra for the calculation of the nuclear 7-ray line emission have been divided into two distinct 
components: 

(1) Standard CRs whose origin and acceleration are widely attributed to massive stars through 
the action of their stellar winds followed by their explosions as supernova. Their composition and 
energy spectra are relatively well known above several GeV per nucleon from direct measurements 
while the extrapolation to lower energies can be done with propagation and demodulation calcula- 
tions. An important indirect constraint on the CR spectra above several hundred MeV per nucleon 
is the diffuse interstellar high-energy 7-ray emission due to vr'' production and decay. This emission 
has been observed frorn sever al Galactic sources by Fermi- LAT, as for ex ample from the local gas 
and dust (jAbdo et al. I ((20091)) and from the inner Galaxy (jStrong I ()201ll )) which are of particular 
interest for the present study. These standard CRs produce an ionization rate well below the mean 
ionization rate deduced from Hj^ column densities observed in diffuse molecular clouds. 

(2) A distinct low-energy component added to the standard CRs to reach the observed ioniza- 
tion rate. For this we adopt the values deduced from recent observations of the molecule Hg" in dif- 
fuse interstellar clouds poi nt ing to a mean ionization rat e of mo l ecular hydrogen of G = 3-4x10"^^ 



s-i (iMcCall et al. I (|2003l ): llndriolo. Fields, fc McCall I (l2009l l: llndriolo fc McCall I (|2012l )l. An- 



other requirement to this component is that the total interstellar CR spectra stay compatible 
with the observed CR spectra above 1 GeV per nucleon ensuring also that the high-energy 7-ray 
emissions stay compatible with the above cited Fermi-LAT observations. 



2.1. Standard CRs 



In the most important energy band for comparison with Fermi-LAT data, between about 1 
GeV and a few hundred GeV per nucleon, data from some ten balloon-borne and satellite-borne 
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instruments exist. We adopted the most recent observations of CR fluxes from the PAMELA 
detector onboard the Russian Resurs-DKl spacecraft w hich provides high-pre cision data for protons 
and hehum nuclei in the rigidity range R ~ 1-1000 GV (jAdriani et al. I (j201ll )). Their proton fluxes 
agree in the kinetic energy range above E 30 GeV, where solar modulation effects should be 
negligible, with other recent data typically at the 20% level or better. The agreement is slightly 
worse for the heliurn fluxes, where differences up to 40% can be observed (see e.g. Fig. 1 in 



Adriani et al. 



(2011 



))• 



To obtain the local interstellar (LIS) proton and helium fluxes the observed fluxes have been 
demodulated with the forc e- field model of iGleeson &: Axford I (|l968l ) with a solar modulation pa- 
rameter of = 450 MV ([Adriani et al. I (|201ll )). This yields LIS proton fluxes in the range E 
= 0.89-1010 GeV and LIS helium fluxes in the range E = 0.36-505 GeV nucleon-^ For the LIS 
helium fluxes we supposed a pure a-particle composition, ignoring an eventual small contribution 
of ^He nuclei. 

The extrapolation of the proton fluxes to lower energies was done with a standard leaky-box 
model for Galactic CR pro pagation. An unbroken p o wer law in parti cle momentum as it results from 
diffuse shock acceleration ( Jones Sz Ellison I ( 199 ll ): iDermer I ( 20121 )) has been used for the source 
spectrum: F(p) oc p~'^. Using dF/dE = dF/dp dp/dE = dF/dp 1/v, the corresponding energy 
spectrum is F(E) = Fq p~'^/f3, where v, p mean the particle velocity and momentum and f3 = 
v/c. Fq and the spectral index s were then adjusted such that the propagated spectrum reproduces 
the LIS proton fluxes obtained from the PAMELA measurements. In the pro pagation calc u lation 



the escape length distribution was taken fror n the disk-halo diffusion m odel of lJones et al. I ()200ll ) 



and the proton inelastic cross sections from iMoskalenko et al. I (|2002l ). The best adjustment was 
achieved with a source index s = 2.35. Figure [1] (left panel) shows the propagated proton spectrum 
together with the Pa, mela da ta and some other rece nt measure ments (IMenn et al. (|200d ) [IMAX]; 
Shikaze et al~l [Bess]: lAlcaraz et al. I (l2000al lbl) [AMS] : boezio et al. I ^199% [Caprice]). 



For the LIS helium spectrum we took the demodulated helium data of PAMELA from E = 
0.36 - 300 GeV nucleon"^ only, the two data points above that energy suffering a large uncertainty. 
Extrapolation of the LIS helium spectra to lower energies was done as in the proton case with the 
leaky-box model with the same para meters for the escape lengths. The inel astic cross section for a 



+ ^He was taken from the model of iTripathi. Cucinotta. fc Wilson I (jl999l ) except below E^ = 15 



MeV where a transition to the experimental values around 10 MeV was operated. Here the best 
adjustment of the Pamela He data in the range E = 1-100 GeV nucleon"^ was found with a source 
index s = 2.32. As shown in Fig. [H the propagated helium spectrum misses the relatively steep 
fall of the LIS spectrum below about 0.4 GeV per nucleon. This part is certainly subject to large 
demodulation uncertainties, the force-field model b eing only an approxim ation and may especially 
fail at these low energies. In the case of electrons, IStrong et al. I ( 201ll ) showed recently through 



constraints on the interstellar electron spectrum from synchrotron radiation that modulation may 
be significantly overestimated for them and that their true interstellar spectrum shows a sharp turn 
over below ~1 GeV. However, as this part has only a minor influence on high-energy 7-ray emission 
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calculations for standard CRs we simply opted for the values of the propagated spectruin . Above 



300 GeV per nucleon we used the data of the ATIC-2 experiment (jPanov et al. I (j2009l )) which 
have smaller uncertainties than the Pamela data there and agree with them at lower energies. The 
spectra of heavier n uclei were taken identical to the helium spectrum with relative CR abundances 
given in table 9.1 of iLongair I (|l992l ). averaged over the three energy bands. 



Wit h these CR spectra the 7-ray emissivity due to tt^ decay has been calculated with the 
model of lKamae et al. I ( 20061 ) for the differential 7-ray emission cross section in the vr'^ — > 27 decay. 



of 


Dermer 


( 


1986). For the 


Murphv. Dermer. & Ramatv 



has been ob tained by multiplying the p + p cross section with the multiplication factor uiHep of 



Mori I ()2009l ) and a linear interpolation was applied in the energy gap. All other reaction channels 
(a + ^He and p,a + HI) were calculated with Mori's multiplication factors relative to the p + 
'^He cross section. Here, HI stands for elements C, N, O, Ne, Mg, Al, Si, S, Ar, Ca, Fe, and Ni. 
For the ambient local interstellar composition, solar photospheric abundances have been used. In 
the CR composition we added ^^B and ^^Cr and ^^Mn as representatives of Sc-Mn, which are 
relatively abundant in cosmic rays. The result of the calcu lation is shown in F ig. [2] together with 



the Fermi-LAT data for the local diffuse 7-ray emissivity of lAbdo et al. I (j2009l ) and their estimate 



of the bremsstrahlung component. The calculated curve reproduces remarkab ly well the observed 



emissivity in shape and absolute normalization, very close to the calculation of lAbdo et al. I (|2009l ) 
albeit with slightly different CR spectra and composition. 

The Galactic CR electron spectrum is constrained by its electromagnetic emissions induced 
during p ropagation in the Galaxy and by direc t data from many balloon and satellite experiments 
(see e.g. IStrong. Moskalenko fc Reimer I (j2009l )). The emissions span a wide range from the Ra- 
dio band for the synchroton radiation to MeV and TeV energies for bremsstrahlung and inverse 
Compton scattering. For the standard CR electron spectrum we rely on results of the GALPROP 
model, which has been developed to calculate numerically the Gala ctic propagation of had r ons an d 
leptons including secondary-particle production and propagation (jStrong fc Moskalenko I (119981 )) . 
The model makes also detailed calculations of the electromagnetic emissions which are induced by 
primary a nd secondary particles. From the different CR electron spectra and propagation models 
studied in IStrong &: Moskalenko I (|2000l ) . we adopt the conventional model. It fits the direct data 
above several GeV and agrees simultaneously with the Galactic synchrotron emission. The latter 
constraint requires a spectral break of the CR elec tron spectrum below several GeV, which was 



confirmed in a recent study by IStrong et al. I (|2011 



2.2. Low-energy part 



The above described standard hadronic CRs yield an ionization rate for molecular hydrogen 
of C2 = 4.3x10^^^ s~^ integrated in the CR energy range E = 0.002 - 10^ GeV nucleon"^, about a 
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fact or of ten below the most rec e nt ob served mean rates in diffuse mo lecular clouds C2 = 4 x 10 



(Indriolo. Fields. & McCall 


( 


2009 


)) and C2 = 3.5+^3 


a direct comparison with 


Indriolo. Fields. & McCah 



ixlO"^^ s"^ 



(jlndriolo fc McCall I ^2012 )) 



energy CR (LECR) component above 2 MeV per nucleon to account for a high (2, we adopted in 
the following for calculations and figures the first value. 



Using the convention al model electron spectrum of 



Strong Moskalenko I (|200Ci ) extrapolated 



down into the eV range, iPadovani. Galli &: Glassgold I (j201ll ) found a negligible contribution of 
electrons to the ionization rate in molecular clouds (C2 ~10~^^ There could still be an 

important CR electron component below 200 MeV to produce the observed ionizatioii rate. Such 
a LECR electron co mponent has for example been studied by IStrong &: Moskalenko I (j200d ) and 
Porter et al. (j2008l ) in order to reproduce the Comptel MeV-excess. But this is out of the scope of 
this paper and we will concentrate in the following on the hadronic component. Supposing that the 
extra ionization in these clouds is entirely due to low-energy nuclei we explored a range of spectra 
and compositions for such a LECR component. 



Scherer et al. I (j2008al ) proposed that a CR population of anomalous cosmic rays (ACRs) accel- 
erated in astrospheres of stars similar to the Sun makes up a significant part of the interstellar CR 
spectrum below E = 0.3 GeV, and when adding the contributions of oth er stars those ACRs may be 



the do minant component of LECRs. Such a scenario was discussed by llndriolo. Fields. &: McCall 



(|2009l ) for their low-energy power-law spectrum (dubbed "carrot" spectrum) with power-law index 
and normalization adjusted to produce the observed ionization rate. We adopted this model in the 
first scenario for our low-energy part of the interstellar spectrum. 

The spectra and composition were derived in the following way: LE CR source spec t ra wer e 



( 2008a 



constructed from the ACR proton spectra at the helio- and astropause of IScherer et al. 
which were constrained by recent data of the Voyager spacecraft. These spectra are well described 
by a power-law function with exponential cutoff at Ec: F(E) = Fq E~^ eT^I^" with a spectral index 
of s = 2.4. We explored ACR spectra for this index and an upper and lower limit of s = 2.7 and s = 
2.0, respectively. For a particles and heavie r nuclei the same index was used and the cutoff energy 
was taken following the species scaling of Isteenbergl (|200nl ): Ec(A,Q)/Ec(p) = (A/Q)-'^^/^^+'^^ 
with energies expressed in energy per nucleon and A, Q being the mass number and the charge 
state, respectively. For sake of simplicity we used 7 = 1 and considered only singly-charged ions, 
which represent the majority of ACR ions in the heliosphere, reducing the scaling to Ec(A,Z)/Ec(p) 
= A~^. Ec was varied from 5 MeV to an upper limit, where the particle fiuxes are still compatible 
with the observed CR helium and proton fiuxes above 0.5 and 1 GeV per nucleon, respectively. The 
upper limits were found to be Ec = 600, 1200 and 4800 MeV for spectral indices s = 2.0, 2.4 and 



2.7, respectively. The ACR composition was taken from lCummings. Stone, fc: Steenberg I (l2002l ). 



For the second scenario we used source spectra resulting from shock acceleration (SA-LECRs) 
as for the standard CRs, but with an energy cutoff Ec: F(E) = Fq R~'/P e"^/^-. R is the 
particle rigidity; spectral index s and E^ expressed in energy per nucleon, were taken the same 
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for all species. For the source spectral index we used s = 2.35, which is the index found for the 
standard CRs as our central value, an d s = 2.0 and s = 2.7 f al ling in the typical r ange of models 



of diffuse shock-acceleration (see e.g. iJones &: Ellison I (jl99ll ): ISchlickeiser I (|2002l )). The energy 
cutoff Ec, which in this scenario may be considered arising from a tuncation in space or time in the 
acceleration process was varied from 5 MeV to an upper limit compatible with the observed CR 
helium and proton fluxes above 0.5 and 1 GeV per nucleon, respectively. It was found to be Ec = 
150 MeV for s = 2.7 and Er = 120 MeV for s = 2.35 and s = 2.0. For the com position o f the s ource 



spectra we took the CR source abundances of iLund I (|l989l ) from table 9.1 in iReames I (jl999l ) 



These spectra were then propagated for both scenarios with the leaky-box model to obtain 
the interstellar LECR spectra. All propagated LECR spectra were normalized to produce an 
ionization rate C2 = 3.6x10"^^ s~^ in diffuse clouds, yielding G = 4xl0~^^ s~^ when added to the 



standard CRs. We calculated the ionization rate as llndriolo. Fields. &: McCall I (|2009|), assuming 
in particular that CRs of kinetic energy below a threshold energy of 2 MeV per nucleon cannot 
penetrate diffuse clouds and thus do not affect the ionization rate there. 

Examples of LECR proton and o-particle spectra are shown in Fig. [3] together with the 
standard CRs alone. The curves with s = 2.35, Ec = 120 MeV and s = 2.4, Ec = 1.2 GeV 
for SA-LECRs and ACR-LECRs, respectively, represent upper limits still compatible with the 
demodulated proton and helium data. The ionization rates of LECRs a dded to the standard CRs in 



dense clouds, taking a low-energy threshold of 10 MeV as proposed in llndriolo. Fields. &: McCall 



(|2009l'l. are displayed in Fig. HI Ne arly all spectra yield rates above the recommended value of 
van der Tak fc van Dishoeck I (I2OOOI ) for molecular cloud cores and with the exception of ACR- 
LECRs with s = 2.7 cross their upper limit for the higher cutoff energies. We note, however, 
that the energy threshold of 10 MeV is an approximation based on stopping range calculations 
and certainly subject to considerable uncertainties. Furthermore, mag netic mirroring can decrea se 
the CR ionization rate in molecular cloud cores by a factor of ~3-4 (jPadovani &: Galli I (1201 ih ). 
which woul d bring all rates if not into the recommended range but at least below the upper limit 
discussed in Ivan der Tak &: van Dishoeck I (I2OOOI ) . 



3. Nuclear 7-ray line emission of the inner Galaxy 



Nuclear line emission was calculated for the same reaction channels as for the calculation of 
vr*^ production except the p + p cha nnel which produces n o 7-ray lines. Cross sections have been 



taken from the recent compilation of lMurphv et al. I (|2009l ). which contains cross section excitation 



functions of about 140 relatively strongly produced 7-ray lines by energetic proton and a-particle 
reactions in astrophysical phenomena like solar flares. The excitation functions extend from reaction 
threshold to typically 1 GeV per nucleon which in most cases is amply sufficient for the present 
study where the nuclear 7-ray line emission is dominated by the low-energy part. Above that range, 
constant cross sections were supposed. 
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For all other 7-ray lines we reli ed entirely on calc ulations with the reac tion code TALYS 



(jKoning. Hilaire. fc Duivestiin I (|2008l )) as it was done in iMurphv et al. I (j2009l ). with some slight 



parameter changes for ^^N, ^'^'^^Ne and ^^Si which gave an improved reproduct ion of new experi- 



menta l 7-ray data for proton and a-particle induced reactions with these nuclei (jBenhabiles et al 
(j201ll )). This code works in the range up to 250 MeV for proton reactions and 62.5 MeV per 



nucleon for a-particle reactions. Above these energies a constant cross section was supposed. Lines 
emitted during the decay of radioactive nuclei were included in these calculations but not the 511- 
keV line and continuum emissions resulting from /3+ decay and subsequent positron annihilation. 
Also emission induced by vr "'"-decay positrons was not included in the calculated spectra. 

As described above, the cosmic-ray spectra were normalized such that an ionization rate C2 = 
4xl0~^^ s~^ was obtained and that the calculated high-energy emission d ue to ir^ decay r eproduced 
the Fermi observations of the diffuse 7-ray emission of the inner Galaxy (jStrong I ()201ll )). For the 
ambient medium composition solar abundances with two times solar metallicity were taken. The 
Fermi-LAT data and calculations for ACR-LECR spectra with s = 2.4 and two extreme cases of 
energy cutoffs Ec = 5 MeV and Ec = 1.2 GeV are shown in Fig. [5l When adding the inverse 
Compton component fro m high-energy leptons and the estimated emission due to unresolved point 
sources as presented in IStrong I (j201ll ) the reproduction is quite good for both cases. The vr" 
production is apparently completely dominated by the standard CRs, even for LECRs with Ec = 
1.2 GeV which add only a small component below E ~ 1 GeV. 

Examples of nuclear 7-ray line emission for three cases of ACR-LECRs with s = 2.4 and the 
case with strongest emission s = 2.0, Ec = 120 MeV of SA-LECRs are shown in Fig. [6l Emission 
due to the standard CRs alone, as described in Sect. 2.1, is also shown for comparison. The total 
7-ray line intensity for the latter is dominated by interactions of the heavy-ion content due to their 
enrichment in standard CRs with respect to the interstellar medium. It accounts for about 70% of 
the total nuclear line emission. Some outstanding narrow lines, produced in proton- and a-particle 
induced reactions are nevertheless present. The emissions of ACR-LECRs with cutoff energies Ec 
= 5 MeV and 25 MeV add practically only a narrow-line component, the quasi-continuum being 
completely dominated by the standard CR component of the CR spectra. In fact, this is mainly 
because heavy ions play a completely negligible role in low-energy cutoff ACR-LECRs due to their 
very low cutoff energies {Ec(A,Z)/Ec(p) = A~^; see Sect. 2.2). Their contribution is less than 10"'^ 
at < 8 MeV, ~6% at Ec = 80 MeV and raises to ~50% at Ec = 1.2 GeV. 

The overall emission is stronger for SA-LECRs because more particles above the reaction 
thresholds are present with respect to ACR-LECRs of the same cutoff energy. ACR-LECR spec- 
tra with their steep rise towards low energies contain many low-energy particles that contribute 
strongly to ionization but are very inefficient for 7-ray line production. This is true for the narrow 
lines induced by light ions and slightly more pronounced in the quasi-continuum, which reflects 
furthermore the importance of energetic heavy ions in SA-LECRs. They contribute between ~50% 
at i?c = 5 MeV and ~60% at Ec = 120 MeV to the nuclear 7-ray emission. 
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The strongest narrow lines in decreasing order are the 4.4-MeV hne (essentially from ^^C and 
^^B, produced in reactions with ^^C and ^^O), the 6.1-MeV line (mainly inelastic scattering off 
^^O and spallation reaction products ^^O and ^^N) and the 1.63-MeV line (mainly from inelastic 
scattering reactions with ^"^Ne). Integrated fluxes of these narrow lines and the integrated fluxes in 
energy ranges 1-3 MeV and 3-8 MeV are shown on Fig. [71 Although the continuum fluxes increase 
steadily with cutoff energy, the narrow-line fluxes tend to level or even slightly decline above Ec 
~50 MeV. The latter can be explained by the cross section excitation functions of the dominant 
line-producing reactions which have their maxima generally well below 50 MeV. An exception are 
6.13-MeV fluxes for SA-LECRs which include an important contribution of ^^0 spallation reactions 
with higher thresholds. 

Narrow-lin e and continuum fluxes have also b een calculated with two spectra and compositions 
as proposed by llndriolo. Fields. &: McCall I ([20091) . The flux values for their " carrot" spectrum are 
within 20% of the values for ACR-LECRs with s = 2.7, Ec = 8 MeV while the values for their 
power-law spectrum agree to better than 15% with the ones for SA-LECRs with s = 2.0, Ec = 
14 MeV. The 4.4-MeV a.nd 6.1-MeV narrow line flu x es are smaller by about a factor of ten than 
the fluxes predicted by [indriolo. Fields. &: McCall I ([20091 ) for the central radian of the Galaxy 
(3 30° < 1 < 30°, |b| < 10°) probab ly mainly because of the method of absolute normalization. 
In [indriolo. Fields. &: McCall I ([20091 ) it was based on the column density of the target material 
in the central radian assumed to be 10^^ cm~^. By normalizing th e nuclear line e mission to the 
vr'^-decay flux from the inner Galaxy as determined by Fermi-LAT ([Strong I ([201 ll )) our 4.4-MeV 
and 6.1-MeV line fluxes are independent of the column density. They depend only on the spectral 
shape and composition of CRs and the metallicity of the inner Galaxy. Supposing that the intensity 
of the CRs which are responsible for the bulk of the tt" production {E ~0.5 - 20 GeV per nucleon) 
is the same throughout the Galaxy, the Fermi-LAT data imply a mean column density of about 
2x10^^ cm~^ in the inner Galaxy (300° < 1 < 60°, |b| < 10°). This explains at least a large part 
of the difference. 



4. Discussion 



Our estimations of nuclear cross section uncertainties can roughly be divided into two cat- 
egories. (1) S trong lines whose cro ss sections excitation functions are explicitely listed in the 



compilation of Murphy et al. (j2009l ). They are at least partly based on experimental data whose 
uncertainties are of the order of 15% or less. For the 1.63, 4.4 and 6.1-MeV lines, the bulk of the 
emission is induced by particles in the energy range where a wealth of experimental data exist, 
and the nuclear uncertainties on the narrow-line fluxes should be less than 20%. (2) The weak-line 



quasi-continuum which is taken from TALYS ([Koning. Hilaire. &: Duivestiin I ([2008[ )) calculations. 
The predictions of this code have been extensively compared to experimental 7-ray c r oss se ctions 



in several studies ( [Benhabiles et al. [ ([201 ih [Murohv et al. I (|2009|) [Tatischeff et al. I ([2006|)) and 



a typical agreement of a factor of two or better for individual lines was found. As this quasi- 
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continuum consists of thousands of individual lines, the overall uncertainty on the integrated cross 
section in the two energy bands 1-3 and 3-8 MeV should be smaller than that. A conservative esti- 
mate would be less than 50%. This weak-line emission component is maximal for the high-energy 
CR component (~70%), while for the LECRs it represents ~50% at the highest-energy cutoffs and 
at = 5 MeV drops to ~20% and ~25% for ACR-LECRs and SA-LECRs, respectively. It results 
in nuclear uncertainties of the 1-3 and 3-8 MeV fluxes of less than ~40%. 

Another uncertainty of the narrow-line fluxes comes from the metallicity of the ambient 
medium. Measurements generally agree for the average metallicity a t the Galactocent r ic dis - 
tance of the Sun to be equal or slight l y bel o w the solar system valu e (jPafion &: Cunha I (|2004l ). 
Lemasle et alTI (|2007l '). IPedicelh et al. I (|2009l ^. iRolleston et al. I (12000)). The radial gradients how- 
ever, especially towards the Galactic center (R = 0), appear less well defined. Using the available 
extracted radial gradients, and assuming they stay constant towards th e center, extrapol a tion o f 
element abund ances to the center yie ld values ranging from [0/H] = (jPaflon &: Cunha I (|2004l )) 
to [Fe/H] = 1 (|Pedicelli et al. I jlQ^d )). where [O/H]=loglO((O/H)i?=o/(O/H)0), and (O/H)©, the 
solar system value. The fact that the inner Galaxy largely dominates the 7-ray emission and pos- 
itive metallicity gradients, however, make a metallicity higher than solar reasonable in this case. 
Our choice of taking [X/H] = 0.3, X = C-Ni for the average metallicity towards the inner Galaxy 
is only a guess with a conservative estimated uncertainty of a factor of two. This uncertainty af- 
fects, however, less the quasi-continuum fluxes which are to a great part due to the CR heavy-ion 
component. Adding that to the nuclear uncertainties a typical overall uncertainty of a factor of 
two is expected for the calculated narrow-line fluxes and between 50% and a factor of two for the 
continuum fluxes. 

The explored parameter space for LECR spectra and composition spans a wide range in spec- 
tral forms and composition, from extremely metal-poor for ACRs with low Ec to the metal-rich 
propagated CR source abundances for SA-LECRs. All spectra are compatible with the high-energy 
diffuse emission of the inner Galaxy observed by Fermi-LAT, the ionization rate of diffuse clouds, 
and recent satellite and balloon data of the CR proton and helium spectra. Estimated ionization 
rates of dense clouds, however, are for all CR spectra, except for the low cutoff-energy ACR-LECRs 
somewhat higher than the recommended value. These estimatio ns are very uncertain, ho wever, as 
they neglect the effects of magnetic flelds on the CR propagation. iPadovani &: Galli I (|201ll ) recently 
found that magnetic mirroring always dominates over magnetic focusing, implying a reduction of 
the CR ionization rate in molecular cloud cores by a factor of ~3-4. 

It would of course be very surprising if the LECR component could be described by a single dis- 
tribution with one definite spectral shape, index, cutoff energy and composition, and be exclusively 
composed of nuclei. However, whatever the nature of LECRs, they can probably be approximated 
by a combination of the spectra and compositions studied here and we believe that the nuclear 
7-ray line emission falls into the range which we have explored and which is presented in Figs. [H 
It is therefore interesting to compare our predictions with 7-ray instrument sensitivities. 
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Integrating over the solid angle of the inner Galaxy (300° < 1 < 60°, |b| < 10°) yields flux 
values (in cm~^ s~^) 0.73 times the values presented in Fig. [71 These values are unfortunately far 
below the narrow-line sensitivities of SPI/INTEGRAL for diffuse emissions: ~5 x 10~^ cm~^ 
for the 4.4- and 6.1-MeV lines with F WHM widths AE ^100 keV an d ~2 x 10"^ cm'^ s'^ for 
the 1.63-MeV line with AE ~20 keV (jTeegarden &: Watanabe I (|200g)). An important fraction of 



the 6.1-MeV line may be ver y narrow, if much of the interstellar oxygen is locked up in dust grains 
(jTatischeff &: Kiener I (j2004l )). It would result in a sensitivity gain of a factor of 5 at maximum, 
but even that is still far from the highest estimated fluxes. In summary, no constraints on LECRs 
can be inferred from SPI/INTEGRAL observations of the inner Galaxy. 

A definite detection of CR-induced nuclear 7-ray line emission or stringent constraints on the 
nature of LECRs in the case of non-detection has probably to wait for the next-generation telescopes 
with an expected increase in sensitivity of typically a factor of 10-100. We therefore compare the 
calculated fluxes to instrument proposals to the last call (2010) of ESA's 2015-2025 Cosmic Vision 
program. Of the three instruments with similar estimated sensi tivities in the .1-10 MeV range 
which have been presented we focus on the CAPSiTT proposal (jLebrun et al. I (I2OIOI )) for which 
some of us participated substantially in the performance estimates. 

The most promising narrow line comparing CAPSiTT sensitivities with the predicted fluxes is 
the 4. 4- MeV line, the sensitivity gain for an eventual very-narrow- line component of the 6.1-MeV 
line is less significant for this instrument with an estimated energy resolution at 6 MeV of ~30 keV. 
Assuming a uniform emission in the defined region, the CAPSiTT 3a 5-y survey line sensitivity for 
the 4.4-MeV line (AE = 100 keV FWHM) is ~7 x 10"^ cm~^ s"^ for such an extended emission, 
using only Compton events in the CAPSiTT instrument. This is at the limit for ACR-LECR 
spectra and would allow detection for SA-LECRs with cutoff energies above 15 MeV. Pair-creation 
events, which have not been considered for the proposal to ESA below = 10 MeV can improve 
significantly the sensitivity for such extended emissions, increasing e.g. the effective detection area 
threefold at ii' = 5 MeV. This would allow a detection for nearly all considered LECR scenarios. 



The broadband 1-3 and 3-8 MeV sensitivities in 5 years amount to ~2 x 10 



—9 

cm s 



and 



~1 x 10~ cm~ s~ , respectively. The flux expected from the standard CRs alone is already 
slightly higher than the sensitivity in the 1-3 MeV band and well above the sensitivity in the 3-8 
MeV band. A detection looks therefore very promising for all LECR scenarios. It may, however, 
be difficult to disentangle from other emissions like inverse Compton scattering of CR electrons 
and unidentified point sources, which probably dominate the diffuse flux of the inner Galaxy in 
the 1-10 MeV range. Analysis and discussion of the diffuse interstellar emission as measured by 



SPI/I NTEGRAL and by COMPTEL/CGRO in this energy range can be found in lBouchet et al. 

boul). 



The nuclear 7-ray line emission has probably a latitude a nd longi t ude p roflle similar to that 
of the high-energy diffuse emission observed by Fermi-LAT (jStrong I (|201ll )). This emission is 
approximately constant along 320° < 1 < 40° and concentrated on the Galactic plane, where 
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about 60% of the emission is in a band with b = ±1.5°. Differences may arise because of the 
metalHcity dependence of narrow-hne fluxes, which is neghgible for vr'^-decay emission, essentiahy 
produced in reactions with ambient H and He. Differences of ionization rates in diffuse interstellar 
clouds, sometimes observed for sight lines with srnall an gular separation indicate a non- uniform 



localized distribution of LECRs (jindriolo fc McCall I (j2012l )). This is not surprising since low-energy 



particles have shorter ranges than the standard CRs resulting certainly also in an inhomogeneous 
distribution of LECRs in the inner Galaxy, concentrated around the acceleration sites. This would 
probably affect only slightly the integrated fluxes since they are tied to the mean ionization rate, 
but the spatial distribution of the nuclear line emission would then be different from the high-energy 
emission from vr^ decay. 

These effects would be favorable for an eventual observation by a Compton-imaging telescope 
like CAPSiTT, whose estimated angular resolution at 5 MeV (0.8°) is about a factor of three better 
than that of COMPTEL/CGRO and SPI/INTEGRAL. The point-source sensitivities of CAPSiTT 
in 5 years amount to ^2 x 10^^ cm~^ s~^ and ~3 x 10~^ cm~^ s~^ for the 4.4-MeV line and the 
3-8 MeV band, respectively. Our calculated fluxes of the 4.4-MeV line from the inner Galaxy are 
typically 10-10^ times higher than that. If e.g. the spatial extension of the inhomogeneities is of the 
order of the angular resolution, a detection gets ver y probable. Assumi n g tha t the large scatter of 



observed ionization rates in diffuse clouds (see e.g. Ilndriolo fc McCall I (|2012l )) reflects the LECR 
density, the 4.4-MeV flux variations may reach a factor of 10. This would be very favorable for 
a detection by a futur instrument with the characteristics of CAPSiTT. Our calculated fluxes in 
the 3-8 MeV band are a factor of 10^-10^ higher than the point-source CAPSiTT sensitivity, and 
a concentration of a few percent of this flux in a point source would be detectable. Following the 
scatter of ionization rates, a factor-of-few variations in the fluxes are probable although a non- 
negligible part of the 3-8 MeV band is produced by standard CRs. That inevitably would be 
detected by a futur instrument with good angular resolution. 



5. Conclusion 

The presented nuclear line and continuum fluxes from the inner Galaxy induced by interac- 
tions of CRs in the interstellar medium were calculated for a large variety of LECR spectra and 
composition and it is probable that the real values fall somewhere inside the obtained ranges. The 
total nuclear 7-ray emission from the inner Galaxy is predicted to be in the range (0.1 - 2.0) x 
10"^ cm-2 s"^ for the 4.4-MeV line, (0.1 - 1.0) x 10"^ cm'^ s"^ for the 6.1-MeV line and (0.3 
- 3.7) X 10~^ cm~^ s~^ for the 1.63-MeV line. For the broad-band emissions we find (0.2 - 1.3) 
and (0.3 - 2.1) x 10~^ cm~^ s~^ for the 1-3 and 3-8 MeV bands, respectively. An eventual obser- 
vation of narrow nuclear 7-ray lines with future 7-ray space telescopes based on actually available 
technology appears to be in reach for a large range of LECR scenarios, although challenging for 
scenarios with low energy cutoffs, where fluxes do not exceed a few 10~^ cm~^ s~^ from the inner 
Galaxy. It may also be promising for the continuum, especially the 3-8 MeV continuum flux which 
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exceeds the estimated future instrument sensitivities by up to a factor of 20. This is especiahy 
true for an instrument with good angular resolution which allows the separation of this emission 
from other diffuse emissions and a very inhomogeneous spatial distribution of LECRs in the inner 
Galaxy could not be missed by such an instrument. 

REFERENCES 

Abdo, A. A. et al. 2009, ApJ, 703, 1249 

Ackermann, M. et al. 2012, ApJ, 750, 3 

Adriani, G. et al. 2011, Science, 332, 69 

Alcaraz, J. et al. 2000, Phys. Lett., B 490, 27 

Alcaraz, J. et al. 2000, Phys. Lett., B 494, 193 

Benhabiles-Mezhoud, H. et al. 2011, Phys. Rev. C, C 83, 024603 

Boezio, M. et al. 1999, ApJ, 518, 457 

Bouchet, L., Strong, A. W., Porter, T. A. & Moskalenko, I. V. 2011, ApJ, 739, 29 

Cummings, A. C., Stone, E. C., & Steenberg, C. D. 2002, ApJ, 578, 194 

Daflon, S., & Cunha, K. 2004, ApJ, 617, 1115 

Dermer, C. D. 1986, A&A, 157, 223 

Dermer, C. D. 2012, Phys. Rev. Lett., 109, 091101 

Gleeson, L. J., & Axford, W. I. 1968, ApJ, 154, 1011 

Indriolo, N., Fields, B. D., & McCall, B. J. 2009, ApJ, 694, 257 

Indriolo, N., & McCah, B. J. 2012, ApJ, 745, 91 

Jones, F. C., & Ellison, D. C. 1991, Space Science Reviews, 58, 259 

Jones, F. C., Lukasiak, A., Ptuskin, V., & Webber, W 2001, ApJ, 547, 264 

Kamae, T., Karlsson, N., Mizuno, T., Abe, T., k Koi, T. 2006, ApJ, 647, 692 

Koning, A. J, Hilaire, S., & Duijvestijn, M. C. "TALYS-1.0", Proceedings of the Int. Conf. on 
Nuclear Data for Science and Technology, April 22-27, 2007, Nice, France, ed. O. Bersillon, 
F. Gunsing, E. Bauge, R. Jacqmin, and S. Leray, EDP Sciences, 2008, 211; see also URL: 
http: / / www.talys.eu/documentation/ 



- 14 - 



Lebrun, F. et al. 2010, in Proceedings of the 8th INTEGRAL Workshop "The Restless Gamma-ray 
Universe", PoS (INTEGRAL 2010) 034 

Lemasle, B. et al. 2007, A&A, 467, 283 

Longair, M. S. 1992, High Energy Astrophysics: Volume 1, Particles, photons and their detection. 
Cambridge University Press 

Lund, N. 1989, in C. J. Waddington (ed.). Cosmic Abundances of Matter, A. I. P. Conf. Proc. 183, 
p. Ill 

McCah, B. J. et al. 2003, Nature, 422, 500 
Menn, W. et al. 2000, ApJ, 533, 281 
Mori, M. 2009, Astropart. Phys, 31, 341 

Moskalenko, I. V., Strong, A. W., Ormes, J. F., & Potgieter, M. S. 2002, ApJ, 565, 280 

Murphy, R. J., Dermer, C. D., & Ramaty, R. 1987, ApJS, 63, 721 

Murphy, R. J., Kozlovsky, B., Kiener, J., & Share, G. H. 2009, ApJS, 183, 142 

Rolleston, W. R. J., Smartt, S. J., Dufton, P. L., & Ryans, R. S. I. 2000, A&A, 363, 537 

Padovani, M., GalU, D. 2011, & Glassgold, A. E., 2009, A&A, 501, 619 

Padovani, M., & Galh, D. 2011, A&A, 530, A109 

Panov, A. D., et al. 2009, Bulletin of the Russian Academy of Sciences: Physics, Vol. 73, 564 
Porter, T. A., Moskalenko, I. V., Strong, A. W. Orlando, E., & Bouchet, L. 2008, ApJ, 682, 400 
Pedicelh, S. et al. 2009, A&A, 504, 81 

Ramaty, R., Kozlovsky, B., & Lingenfelter, R. E. 1979, ApJS, 40, 487 
Reames, D. V. 1999, Space Science Reports, 90, 413 

Scherer, K., Fichtner, H., Ferreira, S. E. S., Biisching, I., & Potgieter, M. S. 2008, ApJ, 680, L105 

Scherer, K., Fichtner, H., Heber, B., Ferreira, S. E. S., & Potgieter, M. S. 2008, Adv. Space Res., 
41, 1171 

Schlickeiser, R. 2002, Cosmic Ray Astrophysics, Astronomy and Astrophysics Library, Springer, 
Berlin, Heidelberg, New York 

Shikaze, Y. et al. 2007, Astropart. Phys., 28, 154 



- 15 - 



Steenberg, C. D. 2000, in AIP Conf. Proc. 528, Acceleration and Transport of Energetic Particles 
Observed in the Heliopshere, Indian Wells, California, 5-8 January 2000, (ACE-2000 Sym- 
posium) ed. M.A. Mewaldt, J.R. Jokipii, M.A. Lee, E. Mobius, & T.H. Zurbuchen, Melville, 
New York 

Strong, A. W, & Moskalenko, I. V. 1998, ApJ, 509, 212 

Strong, A. W, Moskalenko, I. V., & Reimer, O., 2000, ApJ, 537, 763 

Strong, A. W, & Moskalenko, I. V. 2009, Proceedings of the 31*^^ ICRC; available at 
http: / /www. mpe.mpg.de/~aws/propagate.html 

Strong, A. W., 2011, Proc. of the 12th ICATPP Conference on Cosmic Rays for Particle and 
Astroparticle Physics, Villa Olmo, Como, Italy, 7-8 October, 2010, Ed. S. Giani, C. Leroy, 
P. G. Rancoita, World Scientific; and arXiv 1011.1381 (2011) 

Strong, A. W., Orlando, E., & Jaffe, T. R. 2011, A&A, 534, A54 

Tatischeff, V. &: Kiener, J. 2004, New Astronomy Reviews, 48, 99 

Tatischeff, V., Kozlovsky, B., Kiener, J., & Murphy, R. J. 2006, ApJS, 165, 606 

Tatischeff, & Kiener, J. 2011, Mem. Soc. Astron. Ital, 82, 903 

Teegarden, B. J. & Watanabe, K. 2006, ApJ, 646, 965 

Tripathi, R. K., Cucinotta, F. A., & Wilson, J. W. 1999, Nucl. Instr. Meth. Phys. Res. B, 155, 349 
van der Tak F. F. S. & van Dishoeck, E. F. 2000, A&A, 358, L79 



This preprint was prepared with the AAS lATp^X macros v5.2. 



- 16 - 




Fig. 1. — Left: LIS proton fluxes from demodulated data of five recent experiments and propagated 
proton spectrum from the leaky-box model. The modulation potential for the force-field approxi- 
mation used to demodulate the different data sets are indicated in the left figure. Right: LIS helium 
fluxes from the same experiments and from ATIC-2. No demodulation was applied to the ATIC-2 
data because they start above 10 GeV per nucleon, where modulation should be negligible. 
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Fig. 2. — 7-ray emissivity as a function of energy fro m the local atomic h ydrogen gas. Fermi-LAT 
data and the curve for e~ bremsstrahlung are from lAbdo et al. I ( 20091 ). The emission from vr*^ 
decay is calculated with the standard CR composition and spectra as described in Sect. 2.1. 
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Fig. 3. — Interstellar proton (black lines) and a-particle (blue lines) s pectra includin g a low- energy 
component from propagated source spectra of the model for ACRs of lScherer et al. I (j2008al ) (solid 
lines) with s = 2.4, = 1.2 GeV and SA-LECRs (dashed hues) with s = 2.35, = 120 MeV. The 
standard CR spectra from demodulated balloon and satellite data and propagation calculations are 
shown with the dotted lines. Demodulated data from recent observations are represented by the 
symbols (see Fig. [T]and text for more details). 
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Fig. 4. — Calculated ionization rates of cosmic rays in dense molecular clouds supposing that 
particles with energies below 10 MeV per nucleon do not penetrate these places. Red symbols 
(connected by the dashed lines) show the values for SA-LECRs with spectral indices s = 2.0 
(triangles), s = 2.35 (squares) and s = 2.7 (circles), blue symbols (connected by the full lines) the 
values for ACR-LECRs, s = 2.0 (triangles), s = 2.4 (squares) and s = 2.7 (circles). The ionization 
rate of standard CRs (0.35xl0~^^ s~^) is added. The dash ed line and the hatched area show the 
recommended value of Ivan der Tak &: van Dishoeck I (j2000l ) for the cosmic-ray ionization rate and 
its uncertainty in dense molecular cloud cores 

((^^^ = (0.28±0.14) X lO-^'^s"^). The dotted line 

represents their upper limit (~1.3xl0~^^ s~^). 
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Fig. 5. — Fermi-LAT d ata of the diffuse 7-ray emission of the inner Galaxy (300° < 1 < 60°, | b 
< 10°) (jStrong I (|201lh ) are shown by the black triangles. Solid lines: Calculated - r-ray emissions 



from TT ^ decay with CR spectra including a LECR component following the model of lScherer et al. 
(l2008al l for ACRs with s = 2.4, = 1200 MeV (blue) and s = 2.4, Ec = 8 MeV (red). The curve 
for standard CRs is indistinguishable from the red curve, since for Ec = 8 MeV the ir^ production 
is entirely due to the standard CRs. The dashed lines show the total emissi on with contri butions 



of inverse Compton scattering and unresolved point sources as estimated in IStrong 
to the calculated vr'^-decay emission. 



(I2OIII) added 
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Fig. 6. — Calculated nuclear 7-ray line emissions f r om the inner Galaxy for CRs with ACR-LECR 
components following the model of lScherer et al. I ()2008al ) with s = 2.4, Ec = 5, 25 and 1200 MeV 
(magenta, red and green lines, resp.) and SA-LECR with s = 2.0 and = 120 MeV (blue line). 
The emission due to the standard CR component alone is shown by the dashed black line. 
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Fig. 7. — Symbols show the fluxes obtained for ACR-LECRs (left figure) with source spectral 
index of s = 2.4 and s = 2.35 for SA-LECRs (right figure), added to the fluxes of the standard CR 
spectrum: integrated narrow-line fluxes of the 1.63-MeV (filled circles), 6.1-MeV (triangles) and 
4.4-MeV lines (filled squares) and integrated fluxes in the 1-3 MeV range (open circles) and 3-8 
MeV range (open squares) as a function of the cutoff energy of the LECR component. The filled 
and open areas show the flux ranges obtained with spectral indices s = 2.7 (lower flux limit) and s 
= 2.0 (higher flux limit) for the LECR source spectra. The values at Ec = 2.5 MeV, connected by 
dotted lines to the LECR values at 5 MeV correspond to the standard CR component alone. 



